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ABSTRACT 



Fundamental concepts have been formulated for the mechanical behavior 
of isotropic, elastic materials that dewet. This has been accomplished 
through an examination of the stability concepts underlying the classical 
inviscid theory of plasticity. From a phenomenological viewpoint, the 
two theories differ in the nature of unloading, the dilatation behavior, 
the form taken by the dewetting criterion, which is analogous to the 
yield criterion of plasticity theory. Some restrictions to an allowable 
functional form of the dewetting criterion are developed, and a specific 
criterion compatible with these restrictions is suggested. A time- 
temperature superposition method is developed, using the proposed criter- 
ion coupled with reaction rate theory. Stress analysis for materials 
which exhibit anisotropy after dewetting is discussed, and sample 
problem solutions are given for quasi-elastic pressurization and slow 
cool-down of an encased hollow cylinder in plane strain, in which the 
extent of dewetting increases with time. 



This work was partially supported by Naval Weapons Center, China Lake, 
California. 



PREFACE 



This work was supported by the Naval Weapons Center, China Lake, 
California, which support is gratefully acknowledged. It represents a 
continuation of work previously reported in NPS-57L18121A and NPS-57L1981A. 
A portion of the first section of this report makes a correction of one of 
the developments in NPS-57L1981A, concerning normality of the dewetted 
strain increment vector. 

This constitutes a progress report outlining the efforts that have 
produced successful results in the continuing work of building a theory 
for performing fracture analyses on dewet table materials. 
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AN ISOTROPIC THEORY FOR DEWETTABLE SOLIDS 



G. H. Lindsey 
J . E . Wood 

Naval Postgraduate School 



INTRODUCTION 



Composite propellants consist of hard, finely ground solid granules 
embedded in a soft polymeric fuel binder in such a way that each particle 
is completely encased in a polymeric shell. A chemical bond exists 
between the polymeric binder and the filler particles, which may constitute 
as much as 75% to 78% of the volume of the material. 

The composite nature of solid propellants provides a capacity for an 
unusual stress-strain behavior, which is produced by a process of 
separation of the basic elements of the material. As external loads are 
applied to the composite, stresses are produced in the bonds between 
binder and filler . If the stresses are sufficiently large, either the 
adhesive bonds at the interface are broken, or the cohesive bonds near the 
interface are broken, and the binder eventually pulls away from the 
particle leaving a void, or vacuole, surrounding it. This is a dissipative 
process, and it results in a change in the material properties of the 
propellant. The phenomenon is called dewetting, and through the manifested 
change in material properties, it contributes to the nonlinear behavior 
of the material. 

In order to analyze structures fabricated from these materials, it is 
necessary to develop a mathematical theory of material behavior which 
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accounts for this nonlinearity. Several idealizations of dewetting are 
proposed in an effort to gain insight into this type of behavior. A 
phenomenological approach based upon a continuum is pursued here; wherein 
the local details of filler particle shape and size, the dissimilarity 
between binder and filler material, as well as the existence of voids 
partially or totally surrounding the filler particles is disregarded. 

Dewetting is examined in the restricted sense of isothermal, 
quasielastic behavior; i.e. viscoelastic effects are excluded. The stress- 
strain law, both loading and unloading, is then further linearized for 
mathematical simplicity. The simplified behavior is believed to describe 
many of the essential characteristics of dewetting materials, and from it, 
information is obtained concerning the dewetted constitutive law and the 
dewetting criterion. 

To accomplish this, the material is examined in light of the stability 
postulate of Drucker, which has contributed much to plasticity theory. 

From the postulate several interesting variations from plasticity are 
developed for dewetting materials; however similarities with this estab- 
lished theory make comparisons with it inevitable. For this reason 
whenever possible the similarities and differences between the two theories 
are discussed throughout the development. 

A final idealization should be noted. The process of dewetting 
under virtually all applied loads will produce an anisotropic material. 

The approach presented herein assumes isotropy, both before and after 
dewetting, and is intended to be a building block of a more general theory 
to follow. 
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NONLINEAR BEHAVIOR IN UNIAXIAL TENSION 



A specimen of composite propellant subjected to uniaxial tension 
typically possesses a nonlinear stress-strain relationship. It is 
virtually linear up to the onset of dewetting, but at that point the 
reinforcement of the binder around the particles is reduced, and vacuoles 
are formed leading to the manifestation of dilatation in the propellant. 
The next increment of stress produces a correspondingly larger increment 
of strain due to the reduced modulus resulting from the loss of constraint 
through dewetting; thus the curve bends over and becomes nonlinear. 

This process continues as more and more bonds are broken over a greater 
percentage of the surface of each particle, producing an increase in the 
volume change and a decrease in the modulus. 

Idealization of Stress-Strain Law 

From the shape of the curve alone, it is impossible to ascertain 
precisely the nature of the nonlinear phenomenon being observed. For 
instance nonlinear elastic or plastic behavior are indistinguishible 
from dewetting during loading; however, these three are significantly 
different upon examination of the unloading behavior. Representative 
forms of actual laboratory traces are given in Figures 1 and 2, which 
show the real material behavior, including hysteresis. The delineation 
between the two is quite clear, but such complex histories are difficult 
to describe mathematically. It has become customary to idealize Figure 
1 by the representation shown in Figure 3. Although there is no real 
material that behaves exactly as depicted, the elements of actual behavior 
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are included, and the idealization has served a very useful function in 
engineering analysis. It is proposed that dewetting be idealized as 
shown in Figure 4, where unloading is elastic with a smaller modulus 
produced by the breakdown of the material during dewetting. 

Within the framework of this idealization, it can be conceived that 
as a material is loaded in uniaxial tension, it behaves elastically up to 
the onset of dewetting. If at this point the load is held constant, 
energy is dissipated within the system by the creation of new surface at 
the particles. The body, holding the same stress level, now accomodates 
itself to a new strain level which is incrementally larger than before 
dewetting. As the load is increased the process is repeated many times 
as indicated in Figure 5, and in the limit of continuous loading, the 
stepped curve approaches a continuous curve as dewetting occurs gradually 
over a range of stress levels. This gives the nonlinear character to the 
stress-strain curve. 

With reference to Figure 5, below the initial dewetting stress, a° , 
the material behaves elastically with modulus E q . When 0 = 0 ° the 
material is capable of sustaining additional increments of stress with 
a continuous decrease in modulus. With this model, the change in shape 
of the stress-strain curve has the general form of propellant behavior. 

At first the dewetting is slight, and the deviation from E q is minor, but 
as dewetting proceeds, the modulus is continuously reduced. 

When dewetting is completed at o 3 > o°, the modulus is once again 

constant. However, fracture may occur before this happens in many 
materials, particularly in propellants where the flaw filled binder cannot 
withstand such high stress levels and follows a non-linear path to failure. 
At any intermediate point o° < o < o^, unloading is assumed to occur along 
a straight line back to the origin. In substance, loading beyond the 
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linear elastic range has produced a transformation in the microstructure, 
which in turn has produced a new elastic material with a reduced modulus. 
After unloading from the nonlinear range, the stress-strain behavior for 
a second loading is the same as the unloading behavior for the first load 
application. 



DECOMPOSITION OF THE STRAIN 



Because of the nonlinear nature of the material in the dewetted 
region, it is necessary to deal with increments of dewetted strain and 
subsequently sum them by integration. Such a procedure is classified as 
an incremental theory for nonlinear materials, and it should be noted 
that within the assumptions of the quasielastic theory, all nonlinear 
behavior is assumed to be due to dewetting; i.e., all strains are small 
enough to be within the realm of linear deformation theory. 

In the theory of plasticity it is a fundamental assumption, based on 
observation, that the total strain is decomposable into elastic and plastic 
components ; i.e., 

e * e e + (1) 

ij ij ij u; 



For dewet table materials, however, this is not the case. In a sense the 
total strain could be regarded as totally elastic, since it is completely 
recoverable by an elastic law upon unloading, yet inelastic strains do 
occur during loading, which contribute to the total strain. In Figure 6, 
the dewetted increment of strain is decomposed for uniaxial tension. The 
resulting total incremental contributions and the corresponding portions 
of strain are uniquely calculable by integration. 
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To generalize, the conventional assumption of strain decomposition as 



expressed in equation (1) is replaced in the present development by 
decomposition of the incremental strain, 
de . 



, e , , d 

dc.. + de 



while using the definition of inelastic increment of strain employed 
previously by Palmer, Maier and Drucker^ with reference to a model of 
elastic-plastic behavior. As adapted here it states: 

The infinitesimal increment of dewetting strain 
de^ is the change in strain at produced by 
the application and removal of an infinitesimal 
stress increment do^ . 



( 2 ) 



CONSTITUTIVE LAW 



The dewetted constitutive law is obtainable by applying the above 
definition of dewetting strain to the proposed model of material behavior. 

To generalize the one dimensional dewetting behavior to a multiaxial state 
of stress and strain, it is convenient to introduce a space in which the 
coordinates of a point are the components of principal stress, o^. The 
dewetting stress in tension is then a point on a dewetting surface f(o°) * 0, 
which is assumed to exist for the type of materials considered here. 
Generalizing further, there exists about the origin a domain called the 
elastic range which is defined as the totality of elastic states for which 
a change da^ produces no dewetting. If, however, for some do^, it is 
possible to produce a nonelastic response, the material is in a dewetted 
state. The totality of these dewetted states bonds the elastic range 
and defines the dewetting surface. 
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To develop the constitutive law, it is convenient to decompose the 
stress and strain into deviatoric and dilatational components. Any load 



path in principal stress space is also decomposable into a deviatoric 
and a dilatational portion. Consider a point on the dewetting surface at 



stress is related to the strain through the initial elastic constant G o< 



involves both deviatoric and dilatational components. However by the 
isotropic assumpt ion, the stress and strain components can be decomposed 
into dilatational and deviatoric parts after dewetting as well as before. 
By equating (3) and (4) a deviatoric dewetted stress-strain can be for- 
mulated 



this space is along the line of a vector extending from the origin to the 
point of intersection of the load path with the dewetting surface. In 
nine dimensional stress space, the deviatoric constitutive law becomes. 



S°, which is reached with purely deviatoric loading. (See Figure 7) The 




(3) 



A stress cycle, producing dewetting and returning to S°, creates a 
dewetting strain increment, de^, which is related to S° through the new 
modulus G. 



S° = 2G[e° + de^] 



(4) 





(5) 



where 



G = G - dG 



o 



This relationship shows that de*? is in the same direction as S°, which in 
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( 6 ) 



A similar relationship for a purely dilatational loading to the 



dewetting surface at S , yields 



A d _ S_ 

de 3K K 
o 



dK 



( 7 ) 



For a general loading, the dewet ted strain increment is the appropriate sum 
of the two basic parts (6) and (7) 



de d = _iidG + s_cnc 5 

ij 2G G 9K K ij 
J o o 



( 8 ) 



The elastic portion of the strain can also be determined from the 
cycle. Consider the material just after adding the stress dS^ at 
A new elastic material has been generated, and the governing relationship 
between stress and strain at that point is expressed by, 



s i + ds i = 2Gte ° + de i ] 



(9) 



Comparing with equation (4) , the nine dimensional form of the relationship 
is 



dSj 



2Gde* 



'ij ij 

The counterpart expression for dilatation is 

dS - 3Kde e 

The complete dewetted constitutive law in incremental form can be 
written as a sum of (8) , (10) , and (11) 



( 10 ) 



( 11 ) 



de 



ij 






9K ij 2G G 9K K ij 



( 12 ) 



DRUCKER'S HYPOTHESIS 



In 1952 it was shown by Drucker that most of the fundamental 
results of inviscid plasticity theory were derivable from a single assumption 
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regarding the positive definiteness of a certain work-like quantity, 

3 

Drucker extended this concept in 1959 to include inelastic materials with 
rising stress-strain curves, which he termed stable inelastic materials. 
The concept upon which Drucker based his work has come to be referred to 
generally as Drucker* s Postulate or Drucker *s Hypothesis. 

As to the mathematical model of a work-hardening, elastic-plastic 



material, the hypothesis leads to the expression 



<o ij ■ dE ij + 1/2 d °ij dE ij * ° 



(13) 



o * 

where is any state of stress lying on the yield surface, is any 



state of stress inside the current yield surface, and de^ is th e increment 
of plastic strain produced by an outward pointing stress increment do 
at o°_. , subject to the assumptions: 



ij 



1. Continuous variation with plastic strain of shape and position 
of the yield surface. 

2. Continuity and uniqueness (reference) of the stress increment 
do, . for an increment of strain de. 



ij 



'ij* 



DO * 

3. Dependence of de^ upon o^ but independence of o^ 

4. The yield surface is smooth. 

There follows from (13) the proofs of 

(i) convexity of the yield surface 
(ii) Normality of the increment of plastic strain de^ (° r the rate 
of plastic strain e ^j) to the yield surface, which leads 
subsequently to the yielded constitutive law. 

Observing the parallels between yield surface and dewetting surface 
in stress space, between increment of plastic strain de^ and increment 
of dewetting strain de^ > and the assumptions associated with their 
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respective behavior, it is natural to investigate the application of 
Drucker's Postalate to dewetting materials. Whereas the inequality (13) 
follows almost by inspection from Drucker’s stability postulate, the 
corresponding expression for the dewetting model requires additional 
attention due to the fact that the elastic response is altered by the 
inelastic deformation and is in this sense nonlinear. A feature similar 
to this in another model (of plastic behavior) has been shown**' to 
invalidate the proof of convexity and indeed lead to a nonconvex yield 
surface. 

For a dewetting material whose behavior is described by the idealiza- 
tion shown in Figure 4, and whose initial state of stress referred to 
above is zero, the postulate is equivalent to the statement that work 
around a closed cycle is non-zero, or that the work done in dewetting 
cannot be recovered. For the general case, where the cycle begins and 
ends at an arbitrary stress level, the postulate has no obvious thermo- 
dynamic interpretation. It may be summarized as follows: 

Consider a body at rest acted upon by a set of boundary 
conditions and body forces. An external agency, distinct 
from the agency causing the existing state of stress, adds 
slowly a set of surface tractions and body forces. The 
displacements of the body will change as the added force 
is applied and equilibrium maintained. The work done by 
the external agency must be positive or zero over a cycle 
of application and removal. 



Consequences of Drucker's Postulate for a One-Dimension Stress State 

With reference to Figure 6 in a tensile test, suppose after dewetting 

initiates at A and loading to B that the load is partially removed along 

* . 

BC to o . With point C as a starting point and with Drucker s hypothesis 

in mind consider the following stress cycle; 
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X A X 

(i) Loading from the existing state of stress a at C to a > a 

at B , where o° is the state at which dewetting is observed to 
begin. 

(ii) Loading from o° to ■ a° + da at D where da is a positive 

increment of stress producing an increment of strain de of 

which a portion is defined as dewetting strain, de^. 

* 

(iii) Unloading to o at point E. 

It is to be noted here that a stress cycle has been completed in 

* 

that the body has been returned to the same stress o , but the strain 

_* * 

state is different. Dewetting has produced a residual strain A « e - e . 
From Figure 6 the work performed by the external agency in the 

♦ 

application and removal of stress is represented by the area CBDE which is 
obviously positive and in agreement with the postulate. 

By inspection, the enclosed area is 

W - 1/2 (de d + A) (o° - a*) + 1/2 dade d £ 0 (14) 

Since this expression can be shown to be positive, it is to be noted that 
Drucker's Hypothesis is consistent With the uniaxial behavior. It is a 
three dimensional expression of this type that is needed to obtain 
fundamental information about the foundations of dewetting theory. 



Drucker's Postulate Applied to Multiaxial Stress States 

* 

Consider a state of stress in a region R, where dewetting has 
* 

not occurred, 0) • The work done by the external agency in 

loading up to the dewetting threshold, f(o°j) • 0, is given by 



W i-H 



<<i - - c l 3 > dv 



(15) 
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After loading to this point, the material is allowed to undergo an 
increment of devetted strain de^ . During the devetting phase, the work, is 

», ■ J R [<«y - V + 1/2 d »lj dc lj] dv < 16 > 

* 

Upon unloading to the same stress, the strain is different, e ±y 

and the work recovered is 



W 3 - -/JimoJj - Oy + + <v]» (17) 

Combining terms, the work for the cycle is W ■ EW^, and from Drucker's 

Postulate W > 0. At any Intermediate state of unloading, a strain is 

reached which is different from that obtained during loading. This 

* * * 



produces a residual strain denoted by A^ ^ , which at is A^ “ “ e ij* 

It is also noted that the volume is arbitrary, permitting the work expression 
to be written as 



W - 1/2 (oj - °ij) (d £ i j + A ij) " 1/2(e ij 

However, from the dewetting model 

, o . , d . . o 

(e lj + d£ lj’ ‘ A ijkl°kl 



'li )d °lj 1 0 



★ £ 
cij " A ijkl°kl 



de ij A ijkl d0 kl 



Using these relationships in (18) gives 



-o^Hde^ + A^) +do„de^i° 



ij ij 



ij 



ij ij 



(18) 

(19a) 

(19b) 

(19c) 

( 20 ) 
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This limit checks to plasticity theory wherein de^ 
equation (20) reduces to equation (13) « 





and 



SHAPE OF DEWETTING SURFACE 



Shape of Cross-section Normal to Hydrostatic Axis 

o * 

Considering loading histories for which o^j - >> do^ * the 

inequality (20) is written 



<V 0 *j> <de u + V i0 

When cast in terms of deviatoric and dilatational components, (21) 
becomes 



( 21 ) 



(s °j ' s ij )(de ij + V + (s ° - s H* 



* /de + V 



) > o 



( 22 ) 



For a purely deviatoric lcfading path to dewetting, S ■ constant, 

<S U ' V (d ‘l + V - 0 

* * 

Recall that V * e^ - e^ is a .vector in principal strain space lying 

* * * 

in the octahedral plane. Since S. .■ 2G e.. and S . . ■ 2Ge . . then 

ij o ij ij ij 



(23) 



v .fiiri.JLl . lildG 

ij 2 LG G J 2G G 



(24) 



o o 

Substituting (24) and (6) into (23), Drucker's Postulate fcay be expreseed 
in terms of the stress coordinates of the cycle. 



< s ij -^><1 + V 

1 dG 

Since ^ > 0, equation (25) reduces to 



(25) 



,o _o 



* * 



S ij S ij " S ij S ij - ° 



(26) 



13 



In stress space^S^ and are vectors, and each terra in (26) is the dot 

product of the vector with itself. Equation (26) states that the magnitude 

o ^ 

of is always greater than or equal to the magnitude of The only 

way that this can be satisfied is for the contour of the cross-section 

of the dewetting criterion surface in the octahedral plane to be circular. 

Thus the shape of the dewetting criterion must exhibit axial symmetry 

about the hydrostatic axis. 



Shape of Cross-section in Plane of Hydrostatic Axis 

* o 

For a purely dilatational loading path to dewetting, and 

equation (22) yields 

(S° - s*)(de d + V) > 0 (27) 

„ ... , d S° dK 

From (7) de = £ 

o 

o * 

When S and S lie in the same quadrant, 

* * 



, d S 
de - 

o 



S dK 

3K K 
o 



* o 

This relationship does not apply when S and S 
since the bulk modulus after dewetting may take 
tension and compression • 

With this restriction, substitution of (7) 



(s °2 _ s *2\ _dK_ > 

b > 3K K - 



(28) 

lie in different quadrants, 
on different values in 

and (28) into (27) yields 

(29) 



o 

First the sign of dK must be determined. If dK is allowed to take on a 

A A 

negative value, then (29) requires | S° | <, |s |. Since S is any dilational 
stress which does not produce dewetting, including the value zero, 
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this requirement could be met only if S° is identically zero, contrary to 
experimental evidence. Thus for any loading path for which S° and S* are 
in the same quadrant, it is concluded that dK must be positive. Then 
(29) can be satisfied only with the condition that the absolute value of 
the dllatational stress component at dewetting must be greater than the 
absolute value of the dllatational stress component prior to dewetting. 

For any purely dllatational loading path, this restriction means that 
dewetting can be produced only by an increase in the magnitude of the 
dllatational stress. In each individual quadrant in principal stress 
space then the dewetting criterion must be a single-valued function of 
the dllatational component of stress. 

Empirical Consideration 

Analysis of the dewetting process on a mlcrosooplc scale Indicates 
that dewetting occurs as a result of tensile failure of either cohesive 
or adhesive bonds in the neighborhood of the filler-binder Interface. 

In the absence of tensile stresses around the filler particles, dewet- 
ting does not occur. Although it is not feasible to determine the precise 
stress values, it is reasonable to assume that imposition of hydrostatic 
pressure alone creates no tensile stresses in the neighborhood of 
filler particles. Hence, dewetting will not be produced as a result 
of hydrostatic pressure alone. This observation suggests that any 
loading path which starts from an undewetted state and proceeds parallel 
to the hydrostatic axis in the direction of increasing hydrostatic 
pressure cannot intersect the dewetting surface. Consequently, the 
dewetting surface must open continually away from the negative hydrostatic 
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axis, or, in a limiting condition, parallel the axis. Results obtained 
from dilatation measurements during uniaxial tension tests under 
superimposed hydrostatic pressure are in agreement with this hypothesis. 



Functional Representation of the Dewetting Surface 

The restrictions imposed by (26) and (29) , together with the 
empirical consideration discussed above, can be satisfied, by a surface 
with the characteristics shown in Figure 8. Now these characteristics 
of the shape of the surface will be translated into restrictions on 
its functional form. 

In a completely general functional form, the dewetting criterion 
may be given in terms of deviatoric and dilatational stress components 
by 

f(S 1J ,S) - 0 (30) 

An equally general form for isotropic materials, which is useful in the 
present case, may be written in terms of stress invariants: 

f(I r J 2 , J 3 ) - 0 (31) 

where 1^ - is the first principal stress invariant, J 2 ■ 

-(SjS 2 + S 2 S 3 + S^S^) is the second principal deviatoric stress invariant, 
and * S^S 2 S 3 is the third principal deviatoric stress invariant. 

The condition derived previously from Drucker's Postulate that the 
dewetting surface was symmetric about the hydrostatic axis means that the 
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surface may be described by the specification of two coordinates - the 

dilatational stress component, and the radius of the cross-section in the 

octahedral plane, which is given by the octahedoal shearing stress. 

Since the dilatational stress is given by S • 1^/3 ant * t * ie octahedral 

2 2 

shearing stress is given by T oct “ J 2 » t * ie dewetting surface may be 

described independently of This is a significant result. 

fd^) - 0 (32) 



Without more experimental data than are presently available for dewetting 
in multiaxial stress states, a specific form of the function in (32) 
cannot be definitely established. As a first estimation of a specific 
dewetting criterion, Lindsey^ has proposed the simplest form, a linear 
combination. 

1/2 

f(I 1 ,J 2 ) - I 1 + aJ^ -3-0 (33) 



a and 6 are constants to be determined from experiment and J 2 is taken to 
the one-half power for dimensional consistency. In terms of the experiment- 
ally determined stress at dewetting in uniaxial tension o u , and in equal 
triaxial tension a and 3 may be evaluated to give 



3o^ - o 



, 1/2 



j r - 3 \ 



(34) 



u 

The dewetting criterion given by (34) describes a surface in principal 
stress space which is a right circular cone whose axis is the hydrostatic 
axis and whose vertex is located at “ °2 “ °3 ” °t* Th e cone °P ens 
continually away from the negative hydrostatic axis. This criterion 
satisfies all the restrictions derived and hypothesized above. At this 
time no claim can be made that (34) is indeed the correct dewetting 
criterion; however it can represent real behavior sufficiently well to be 
useful in analysis. Specific applications are developed in the following 
sections . 
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TIME AND TEMPERATURE EFFECTS 



In the elastic concept of dewetting behavior discussed in the preced- 
ing section, dewetting is considered to occur instantaneously when a 
critical value of stress is reached. With this concept, if the material 
is loaded to produce a stress state less than the critical value, 
dewetting should never occur. In actual experience, however, dewetting 
does not occur instantaneously upon the application of a suitable load. 

The time at which dewetting occurs is found to be dependent upon loading 
rate and temperature, reflecting, in part, viscoelastic material response. 
However, experimental evidence indicates that temperature and loading rate 
influence the time to dewetting even for loading conditions normally 
considered within the range of quasi-elastic behavior. The data reported 
by Fishman and Rinde^ for uniaxial tension tests at strain rates less 

than 1% per minute show such an influence, which incidentally is not 

9 

correlated by use of the WLF equation normally used for time-temperature 
superposition of viscoelastic material response characteristics. 

For further generalization of the concept of a criterion for dewetting, 
the influence of time and temperature on the dewetting process will be 
considered, while retaining the assumption of quasi-elastic stress-strain 
response. The aim of this effort will be to represent the effects which 
are peculiar to the dewetting process separately from viscoelastic effects 
which may also be present in a real material. 

In considering the dewetting surface as a function of time and 
temperature as well as stress, it will first be assumed that the functional 
form does not change with the introduction of these new variables. The 
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dewetting surface developed earlier may then be generalized to represent 
the locus of all stress states for which dewetting will occur at the 
same time, or at the same loading rate, and at the same temperature. 
Lindsey's stress criterion for dewetting, eq. (34), may be rearranged in 
the form 



This has the form in principal stress space of a cone centered on the 
hydrostatic axis with vertex at “ °2 “ °3 “ ° t > the dewetting stress 
in hydrostatic tension. If dewetting is observed for hydrostatic tensile 
stress o at a particular temperature and loading rate or time after 
imposition of load, then the generalized dewetting criterion predicts 
that dewetting will occur at that same time for any other state of stress 
represented on the cone with vertex at o t> when the temperature and type 
of loading are the same. If o is a function of time or temperature, the 
resulting dewetting surface will be displaced along the hydrostatic axis. 
For example, with the vertex at the origin (zero hydrostatic tension), 
the time until the occurrence of dewetting is infinite. As the hydrostatic 
tension is increased, the vertex of the cone moves in the positive 
direction along the hydrostatic axis, so that instantaneous dewetting 
would require an infinite hydrostatic tension. Variation with time or 
temperature of the ratio on the other hand, would change the vertex 

angle of the conical dewetting surface. 

In order to obtain an analytical expression relating time, temper- 
ature, and stress in the dewetting process, the nature of the governing 
molecular phenomena should be considered. For the purpose of modelling 




I 



(34a) 



u 
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response of high polymers, the molecular processes involved in deformation 
under stress are usually divided into three classes-: 

a* Slippage of chain segments which are relatively free to move 
about within the polymer network - a purely viscous response. 

b. Extension of coiled and entangled chains - a viscoelastic 
response which may be modeled by multiple Voigt-Kelvin models 
in series. 

c. Rupture of chains through breakage of chemical bonds. 

Tobolsky and Eyring^ presented a mathematical formulation of these 
responses which they applied to a number of types of deformation processes 
such as creep, stress relaxation, vibration and rupture in polymeric 
materials. They state that the viscosity associated with the sloppage of 
chain segments (a) is much less than the resistance to deformation of 
processes (b) and (c) . They also find that the stresses required for 
bond breakage are significantly higher than for chain extension and thus 
assume that for processes of creep and relaxation, except possibly at 
large deformations or long times, chain extension is the significant 

g 

molecular process occurring. Bueche has shown that consideration of the 

response due to chain elongation for different times and temperatures leads 

9 

to the WLF equation which has been found to represent time-temperature 
superposition for creep and relaxation data in many polymer systems. 

Tobolsky and Eyring published an equation giving the rate of 
breaking bonds in a polymeric thread under any tensile loading: 



I dN 
N dt 




f * applied tensile stress 

\ * distance between equilibrium positions 

N = number of unbroken bonds per unit cross-section area 



( 35 ) 
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AF - apparent energy of activation for the bond breaking process 

K, h, R * Boltzman's constant, Planck's constant, and the universal 
gas constant, respectively. 

When the argument of the hyperbolic sine is sufficiently large, as 
it is under applied loads of the order of 10 psi, and stress and tempera- 
ture are held constant, (35) may be Integrated to obtain the time for 
which the number of remaining bonds per unit area goes from the initial 
value N q , to zero. This is interpreted as the time to failure. 

2N h Xf 

C f " f X eXp ^RT " 2KTN ) 
o o 

Validity of equations (35) and (36) for the prediction of bond 
failures is not restricted to polymeric materials. Equation (36) has been 
found to predict accurately the failure times for uniaxial tensile tests 
of a great variety of materials. Zhurkov 3 ^, for example, found that an 
expression of this form predicted fracture behavior for fifty different 
materials. Including metals and alloys, non-metallic crystals and polymers. 
In an extension to previous work, Graham and Robinson 3 ' 3 ' found that failure 
times for specimens under a time-varying loading history correlated well 
with predictions obtained by numerical Integration of equation (35) for 
such loading paths as constant strain rate, bl-llnear strain rate, and 
multiple cycles. 

The exact molecular behavior involved in the dewetting process is 
not completely clear. Dewetting may occur in either of two modes: an 

adhesive failure, in which the bond between filler and binder is broken; 
or cohesive failure, in which microscopic vacuoles form completely within 
the binder phase and coelesce to form larger vacuoles at the surface of 
the filler particles. Mechanisms of viscous flow may contribute to the 



21 



formation of the microscopic vacuoles, but the coelescence of vacuoles 
and the adhesive failure seem to depend significantly on bond breakage. 

To the extent that this is true, an equation predicting the rate of bond 
breakage such as the Tobolsky-Eyring equation should be useful for predict- 
ing dewetting. One successful application of the Tobolsky-Eyring equation 
to the prediction of time of dewetting initiation was reported by Graham 
and Robinson^ for a polyvinyl chloride plastisol propellant formulation. 

The analysis used by Graham and Robinson for their PVC data considered 
the existence of only one type of failure reaction. They made the further 
simplifying assumption that the ratio A/N^ which appears in the reaction 
rate equation is a linear function of absolute temperature, an assumption 
on which opinion of various investigators is divided. In the relatively 
restricted range of temperatures and stresses in which their data were 
obtained, their representation of the reaction rate theory led to good 
correlation. However, their correlation was not as close for failure times 
in a PBAN propellant formulation, especially at low temperatures and high 
stress rates. 

In order to apply reaction rate theory to the dewetting process 
over the entire range of temperatures and loading rates which may be of 
interest, it may be necessary to consider more than one type of molecular 
reaction. For example, it is possible that some propellant compositions 
may initiate, -dewet ting by a cohesive failure mechanism at high temperatures 
and an adhesive failure mechanism at lower temperatures or higher loading 
rates, with different values of reaction rate parameters associated with 
each process. Also, for loading conditions out of the range of quasi-elastic 
material behavior, the stress in the reaction rate equation should 
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reflect viscoelastic response* One method of including visco-elastic 

response at low temperatures has been proposed by Graham, Henderson and 
12 

Robinson . Finally, temperature dependence of the reaction rate 
parameters should also be experimentally determined. The writer has 



of the absolute temperature yielded a slightly improved fit of the Graham - 

Robinson PBAN failure data and correlation with failure data for a PBAA 

13 

formulation reported by Leeming . In any case, the successful correlations 
obtained by Graham and Robinson indicate that reaction rate theory can 
well represent the relationship between time, temperature, and stress in 
the dewetting process. 

In order to apply reaction rate theory to the multi-axial dewetting 
criterion concept, it is necessary to develop a reaction rate equation 
which reflects the effect of the multi-axial stress state on the dewetting 
process. This may be done in the following manner. 



which dewetting occurs at the same time and temperature in uniaxial tension 

and in hydrostatic tension. Consequently, the dewetting criterion will 

represent the locus of all stress states which will dewet at the same 

time for which a and o were obtained, when the material is at the same 
t u 

temperature. It is reasonable to suppose that the same molecular reactions, 
namely bond breakage , govern dewetting in multiaxial stress states as in 
the uniaxial case. Hence, reaction rate theory would predict that the 
ratio 0/0 would be independent of rate or temperature, although each 
individual term would not. This may be checked experimentally. 



determined that an assumed dependence of the ratio A/N o on the square 




♦q ^ - h - 




(34b) 



where T is the temperature and t^ the time at which dewetting occurs. 

In this expression and will be chosen to represent the stresses for 
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Let the left side of equation (34b) define an equivalent uniaxial 



tensile stress o , This equivalent stress represents the influence of 

e q 

deviatoric and dilatational stress components in the proportions required 
by the dewetting criterion. In the uniaxial tension, case, of course, at 
the onset of dewetting, ® °u^ ,t d^ an< * °u ^ e ^ ne< * implicitly by the 
reaction rate equation (35) . When, for stress states other than uniaxial 
tension, the equivalent uniaxial stress 




is substituted in the reaction rate equation, the shape of the dewetting 
criterion surface will be preserved, as discussed earlier. With this 
substitution, the reaction rate equation can reflect the effects of time 
and temperature on the process of dewetting in multiaxial states of stress. 

If the Tobolsky-Eyring reaction rate equation is chosen as the 
governing equation, then, the criterion for dewetting in multiaxial 
stress states as a function of temperature and time to dewetting is obtained 
by integration of the equation 



a X 



1 dN ZKT f-L . u ( eo_\ \ 

• n dJ ' IT e * p (“Rf) (359) 

to obtain the time required for the number of unbroken bonds to go from 
the initial value N q to zero. In the present development for quasi-elastic 
material behavior, o may be obtained from an elastic analysis. If the 
functional form of the stress dependence of the dewetting criterion is the 
same when viscoelastic response is allowed, then the dewetting criterion 
could be extended to viscoelastic materials by substitution of the equivalent 
stress obtained from a viscoelastic solution. When anisotropic response 
of the dewetted material is allowed, the above criterion must be considered 
valid only for the initial dewetting. 
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Movement of Dewetting Surface 



For a dewettable material which is in a constant state of stress 
at a constant temperature, the devetting criterion in principal stress 
space is then described by (37) . When a particular time to dewetting and 
temperature are specified, is the value defined implicitly by (36). 

The locus of stress states for which o * constant generates a cone 
about the hydrostatic axis which is also the locus of stress states with 
the same time to dewetting. Thus the time-dependent dewetting criterion 
in principal stress space may be described in terms of a moving conical 
surface. At time zero the cone is located at infinity on the hydrostatic 
axis. With increasing time, the cone moves along the hydrostatec axis 
toward the origin. Dewetting is predicted to occur when the existing 
stress state intersects the cone surface, shown schematically in figure 10. 

When stress and temperature are allowed to vary with time, the time 
to dewetting must be calculated by integration of (35) . Thus the rate of 
motion of the dewetting surface is dependent upon the loading path. 



CONSIDERATION OF DEWETTING-INDUCED ANISOTROPY 



The development in the preceding section has considered the case of 
a material which is isotropic and linearly elastic both before and after 
dewetting. As a first step toward greater generality, the affects of 
anisotropy in the dewetted phase of the material will be considered, while 
still restricting response characterization to the elastic or quasi-elastic 
case. 

The initial dewetting criterion is assumed to depend only upon stress, 
and the material is considered to remain isotropic until dewetting occurs. 
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Thus the stress criterion for initiation of dewetting in a previously 

undamaged material may still be given by a function of the form of eq. 

(34) . However, the criterion for initial dewetting will be affected to 

some degree by the type of material response which follows dewetting. 

14 

As White and Drucker have discussed in connection with plasticity theory, 
when anisotropy is introduced as a result of deformation, the subsequent 
loading function (in this case, the dewetting criterion) must depend 
explicitly upon the components of dewetted strain, as well as stress. 
Consequently, the dewetting function is no longer a path-independent function, 
but depends upon the prior history of dewetting in the material. Similarly, 
the dewetted constitutive law is also history-dependent. This may be 
demonstrated with a simple example. A sheet of dewettable material is 
stretched in uniaxial tension until dewetting occurs (Fig 11a) . The onset 
of dewetting is predicted by the previously developed criterion. As numer- 
ous visualization studies have demonstrated, dewetting results in separa- 
tion of binder and filler particles in the direction of strain and produces 
an associated reduction of modulus in that direction. Now if the same 
sheet is next stretched in a different direction (Fig lib) , its response 
will depend upon the orientation with respect to the previous direction 
of dewetting. This new response must be considered in order to predict 
the point at which further dewetting will occur. 

Material Orthotropy 

Some information about the nature of anisotropic material response 
resulting from initial dewetting can be ascertained from a particular 
model. The composite material is represented as an orderly array of 
uniform spherical filler particles within a uniform strain field each 
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surrounded by binder. Figure 12 illustrates a cubical building block of 
such an array, oriented with sides parallel to the principal strain 
directions, with dewetting in the direction of e^. The shaded areas of 
the sphere indicate regions in which the filler-binder bond has been 
broken. Since the faces of the cube are planes of symmetry of the array, 
they remain plane during deformation with their relative motion given by 
the values of the principal strains. 

As has been observed by Elliott^ the areas of dewetted surface 
(shaded areas) and the shapes of the resulting vacuoles are influenced 
by the relative magnitudes of the three strains. Indeed, if the dewetted 
area were specified, an elastic solution could be obtained relating the 
average stresses on the walls of the deformed cube to the strains at the 
walls, giving, in effect, the anisotropic constitutive law of the material. 
If the extent and shape of the dewetted area of the spherical inclusion 
were a known function of the strains which exist when the dewetting 
criterion is reached, then the dewetted constitutive law for the array 
could be calculated from a knowledge of the pre-dewetted strain field. 

At present such a solution is not available. However, if the 
assumption is made that the shape and extent of the dewetted area is 
uniquely determined by the strain field which exists at the initiation 
of dewetting, some information about the resulting anisotropy may be 
deduced. With this assumption, the dewetted material will exhibit a 
geometric anisotropy possessing axes of symmetry paralleling the directions 
of principal strains which exist prior to dewetting. Since, in theory at 
least, the dewetted constitutive law is uniquely determined once the geo- 
metry is specified, the material anisotropy will also have axes of symmetry 
corresponding to the directions of principal strain. If a general 



27 



curvilinear coordinate system were chosen whose coordinate directions 
correspond everywhere with principal strain directions, then those coor- 
dinate directions would be axes of symmetry of the anisotropic dewetted 
material . 

The particular type of anisotropy which exhibits symmetry about three 
orthogonal directions is classified as orthotropy. Thus, with reference 
to a coordinate system corresponding with principal strain directions 
of the pre-dewetted material, a material after initial dewetting may be 
classified as orthotropic. In such a coordinate system, the number of 
independent material moduli required to describe the dewettfcd material is 
reduced from thirty-six to nine. 

When the strain field is uniform within a pre-dewetted structure, 
then the resulting anisotropy will be homogeneous. When a nonuniform 
strain field exists prior to dewetting, then the dewetted material will be 
non-homogeneous with a spatial variation in properties directly related 
to the variation in principal strains. For loading geometries which 
produce complex patterns of principal strains in the isotropic elastic 
solution, analytical evaluation of the stress-strain solution including 
dewetted regions will be complicated. However, one problem of engineering 
interest which permits further simplification and analytical solution is 
a circular cylinder under axisymmetric loading in plane stress or plane 
strain conditions. For this particular geometry the directions of 
principal strain and stress in the isotropic elastic material are parallel 
with the coordinate directions of a standard cylindrical coordinate system. 
Consequently, the ensuing dewetted material will be orthotropic with 
respect to the same coordinate system. Cylindrical orthotropy when reduced 
to the case of plane stress requires only three independent moduli to 
characterize response. The constitutive law is given in terms of 
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engineering moduli by 



< 




and symmetry of the coefficient matrix requires v „/E - v /E . For 

r6 6 8r r 

the case of plane strain, the constitutive law becomes 




1/E r (l - v zr vrz) 




e 




e 



with a corresponding symmetry requirement. Since there is no angular 
variation in strain for the case of axial symmetry in load and geometry, 
the anisotropic moduli will at most be functions of radial position, 
reflecting the radial variation of principal strains. 

Bilinear Approximation 

In the general case, the constitutive laws of a composite material 
both before and after dewetting may be nonlinear functions of strain, 
reflecting material nonlinearity. For many compositions when loaded 
quasi-elastically , the uniaxial stress-strain curves may be closely 
approximated by a bilinear constitutive law, shown schematically in 
Figure 13. For loading the material is considered to have a constant 
modulus prior to dewetting and to transition abruptly to a new constant 
modulus reflecting dewetted behavior when a critical value of strain or 
stress is reached. A model developed by Farris^ relating dilatation to 
the dewetting process shows that this bilinear idealization is equivalent 
to the assumption that all available filler particles dewet simultaneously 
when the critical value of strain is reached. In order to simplify the 
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the following analysis, this bilinear idealization of loading will be 
employed. Of course unloading will be handled in the manner described 
previously. Thus, material tangent moduli will be considered to have 
constant values, both before and after dewetting. 

Example Solutions 

Two particular problems of engineering interest in the solid rocket 

field are internal pressurization and slow cooling of the hollow cylinder 

of composite solid propellant bonded in a flexible case of higher modulus. 

Conditions of plane strain are assumed. Linear isotropic elastic solutions 

for each of these examples are obtained by classical means. A convenient 

tabulation is given in Reference 17. The dewetting criterion, equation 

(34a) , with the dewetting stress in equal triaxial tension o taken as 

twice the dewetting stress in uniaxial tension in agreement with 

Lindsey"*, predicts that dewetting occurs first at the inner bore of the 

solid propellant for both cases under consideration. Thus, the problem to 

be solved is a three-element composite circular cylinder, illustrated 

in Figure 14. The inner ring is composed of dewetted material and is cylin- 

drically orthotropic. The center ring is undewetted propellant and is 

isotropic. The outer layer is the rocket case, and may be either anisotropic 

or isotropic but will be considered isotropic in the present analysis. 

18 

Solutions have been given by Lekhnitskii for a number of problems 

involving cylindrically orthotropic structures loaded in plane stress, 

including one analysis in which material moduli were allowed to vary 

radially according to a specific power law. In addition, a solution 

19 

has been published by Freudenthal for stresses due to slow cooling of an 
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orthotropic cylinder bonded to a rigid case. Using the methods of 
Lekhnitskii, solutions have been obtained for the geometry shown in 
Figure 14 for stresses within the composite cylinder as a function of the 
position of the interface between dewetted and undewetted propellant, for 
conditions of plane stress and plane strain. The effects of radial 
variation of the anisotropic moduli are disregarded. Also, in the absense 
of any information to the contrary, the coefficient of thermal expansion 
for the composite propellant is assumed to be unaltered by the dewetting 
process. The resulting formulas are given in Appendix I. These formulas 
were evaluated with the use of a digital computer for a particular geometry 
and physical properties estimated to be representative of a typical 
steel-encased analog rocket motor. The particular values used and results 
obtained are discussed in Appendix II. 

Using the reaction rate equation of Tobolsky and Eyring with the 
value of equivalent stress defined in (37) , the position of the interface 
between dewetted and undewetted material as a function of time is calculated 
for the previously discussed problems of internal pressurization and slow 
cool-down of an encased cylinder in plane strain. For the problem solution, 
it is assumed that the differential form of the reaction rate equation (35) , 
could be applied for time-dependent values of both stress and temperature. 
The solution proceeds in the following steps: 

1. The isotropic linear elastic solution is examined to determine 
the position at which the highest value of equivalent stress is 
located. In all cases this is the inner boundary. 

2. The reaction rate equation is integrated to determine the time 
at which material at one radial increment from the inner boundary 
dewets. In the case of internal pressurization, a constant stress 
level is used. For the cool-down problem, a constant rate of 
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temperature decrease is considered, and the associated stress 
history is used. 

3. After dewetting takes place at the inner boundary, a new 
stress solution is calculated for a three-element cylinder with the 
inner anisotropic ring one radial increment thick. 

A. Next, the time to dewetting is calculated for a radial position 
two increments out from the inner boundary. For the case of 
internal pressurization, the stress history at this point is: 

a. constant stress as given by the isotropic solution from 
time zero until the calculated time at which the first increment 
dewets . 

b. a second value of constant stress, as given from the three- 
element anisotropic solution. A similar step change in stress 
history is calculated for the cool-down problem. 

5. The point of calculation again moves by one radial increment and 
the time at which the dewetting interface reaches that point is 
calculated as in step A. 

6. The process is continued until the entire cylinder dewets or 
until the equivalent stress in the undewetted area reaches zero. 
Calculations are performed with the aid of a digital computer program. 

No computing problems have been encountered to date. The longest computer 
running time experienced so far have been under one minute. A discussion 
of the specific data used in calculation and the results obtained is 
listed in Appendix III. 
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SUMMARY 



The present development constitutes a first step toward a theory 

of structural analysis of dewettable materials loaded to produce multi- 

axial states of stress. Idealization of the material as isotropic and 

elastic, both before and after dewetting, permits consideration of the 

nonlinear variation of material properties associated with dewetting 

independently from other nonlinear effects associated with tost dewettable 

materials. Treatment, analogous in many respects to plasticity theory, 

yields information about a stress criterion for the initiation of dewetting. 

It is found that the criterion for dewetting is independent of the third 

deviator ic stress invariant, and thus may be described in terms of two 

independent variables: the nean hydrostatic stress, and the octahedral 

shear stress; or, alternatively, the invariants 1^ and J 2 . For illustra- 

1/2 

tive purposes, the criterion al^ + - 8 ■ 0 satisfies the constraints 

of the development, and is in agreement with the limited amount of 
available information concerning dewetting of real materials under quasi- 
elastic loading. 

Effects of time and temperature on the process of dewetting in an 
otherwise elastic material can be considered with the use of a reaction 
rate theory for molecular bond breakage. The dewetting criterion provides 
an equivalent bond stress for generalization of the essentially uniaxial 
reaction rate theory to multiaxial stress states. 

A method of analysis which includes the effects of anisotropy 
induced by dewetting according to the above criterion is developed. It 
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is concluded that upon first dewetting, a material assumes the characteris- 
tics of orthotropy with reference to a general curvilinear coordinate system 
which coincides with the directions of the principal stresses in the pre- 
dewet ted state. Stress solutions in closed firm are given for the specific 
problem of an encased, hollow, circular cylinder of dewettable material 
loaded in plane strain by either internal pressurization or slow cooling. 

With these stress solutions, the time histroy of dewetting within the 
cylinder is calculated by integration of a reaction rate equation. 

Experimental investigation is necessary to justify the choice of 
any specific form of dewetting criterion, to determine the values of the 
parameters used in the reaction rate equation and to evaluate the anisotropic 
moduli of the dewetted material. When these values are obtained, the 
method of analysis discussed herein should provide meaningful quantitative 
results within the limitations imposed by the quasi-elastic idealization. 
Alternatively, qualitative bounds on structural response may be obtained 
by parametric studies. 

The approach used in the development of an analysis for the idealized 
material considered in this report may serve as a basis for a more general 
analysis which should consider effects of viscoelasticity and anisotropy 
in greater detail. 
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FIGURE I . PLASTICITY WITH 
HYSTERESIS 



FIGURE 2. DEWETTING WITH 
HYSTERESIS 





FIGURE 3. PLASTICITY MODEL FIGURE 4. DEWETTING MODEL 
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FIGURE 5. IDEALIZATION OF QUASIELASTIC DEWETTING 
BEHAVIOR 



38 



cr 




€ 



FIGURE 6. STRESS CYCLE FOR A DEWETTING MATERIAL 
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FIGURE 7. PURELY DEVIATORIC LOADING PATH 
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FIGURE 8. ALLOWABLE DEWETTING CRITERION IN 
PRINCIPAL STRESS SFACE 
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FIGURE 10. GENERALIZED DEWETTING QRITERION 
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FIGURE 12. ELEMENT OF DEWETTED MATERIAL MODEL 
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FIGURE 13. BILINEAR IDEALIZATION 




FIGURE 14. ROCKET MOTOR WITH DEWETTING 
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FIGURE 15. RATIO OF STRESS TO INTERNAL PRESSURE 

VS 

RADIAL POSITION 
(ISOTROPIC SOLUTION) 
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FIGURE 16. RATIO OF STRESS TO INTERNAL PRESSURE 

VS 

POSITION OF INTERFACE 
(ANISOTROPIC SOLUTION) 
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FIGURE 17. STRESS VS. RADIAL POSITION FOR 1° F 

TEMPERATURE DROP (ISOTROPIC SOLUTION) 
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FIGURE 18. STRESS VS. INTERFACE POSITION FOR l°F 

TEMPERATURE DROP (ANISOTROPIC SOLUTION) 
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19. TIME TO DEWETTING VS. RADIAL POSITION 
(INTERNAL PRESSURIZATION) 
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FIGURE 20. TIME TO DEWETTING VS. RADIAL POSITION 
(EQUILIBRIUM COOLING) 




FIGURE 21. COOLING SCHEDULE FOR EQUILIBRIUM COOLING 
PROBLEM 
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APPENDIX I 



STRESS FORMULAS 



A. Internally Pressurized Cylinder - Plane Stress 
1. Isotropic solution - no dewetting 





P = internal pressure at a Q . 



q 2 = pressure at interface 
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2. Solution with dewetting. Region 1 is anisotropic 
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K = 

r 

and are the pressures acting at the interfaces a^ and a ^ respectively, 

B. Isothermal Cooling - Plane Stress 

1. Isotropic solution - no dewetting 
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a^, a ^ are the linear coefficients of thermal expansion in regions two and 

three respectively. 

T q is the stress-free temperature. 

2. Anisotropic solution - region 1 has dewetted. 

It is assumed that the coefficient of thermal expansion does not 

change with dewetting. 
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values with the primed quantities listed below. Superscripts refer to 
regions of cylinder. 
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APPENDIX II 



EVALUATION OF STRESS FORMULAS 



Estimation of Material Properties 

For sample calculations of stresses within a hollow encased cylinder, 
the following representative values of material properties were selected. 
For the undewetted isotropic propellant: E 2 " 300 psi 

v 2 - 0.5 

a 2 - 2.21 x 10~ 5 in./in./°F 

For the case (steel) : c 3 x 10^ psi 




c * 3 = 6.5 X 10' 6 in./in./°F 

In order to obtain estimates of material properties of the propellant 
after dewetting, the effect of vacuole formation is considered. Since 
during dewetting the vacuoles form and extend in the direction of the 
maximum tensile stress, which in the cases under consideration is the 
tangential direction, loss of reinforcement is most severe in this direction. 
In radial and tangential directions, the binder-filler bond is still intact, 
consequently modulus values for these directions will be less affected by 
dewetting. As a first approximation it is assumed that the modulus values 
and will not change during dewetting , and will remain equal to the 
isotropic modulus With this assumption the Poisson T s ratios and 

will also remain at their isotropic value 



57 



To estimate the extent to which tangential modulus is changed with 

dewetting, uniaxial tensile test results are considered. Francis and 
20 

Carlton report for a composite propellant investigated by them, an 

approximately 50% reduction in uniaxial Poisson's ratio was associated 

with the dewetted state. For the purposes of the present calculation, the 

same reduction is assumed for the corresponding ratio * v 0 / 2* 

dr z 

Similarly v = ^^th these values determined, the rest of the 

desired moduli may be obtained from requirements of symmetry in the modulus 
matrix. In summary, the following values are chosen for anisotropic 
properties of the dewetted material. 
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For the plane strain solution (see appendix I) the moduli are replaced 
with primed values given by 
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In the absence of any information to the contrary, it is assumed 
that the coefficient of thermal expansion remains isotropic and unchanged 
from the value of the undewetted material* 

Geometry 

A particular set of radius ratios was chosen corresponding to the 
dimensions of an analogue motor for which some test data are available. 
These ratios are 
a 

— - 0.2465 
a 2 

a 2 

— = 0.934 
a 3 



Results 

A computer program was written to calculate the stresses developed at 
various points within the propellant as a function of the position of the 
interface between anisotropic and isotropic propellant per pound of internal 
pressure. Results are shown in Figures 15 and 16. With the high modulus 
steel liner, the isotropic solution is very nearly uniform hydrostatic 
pressure, reflecting the choice of an incompressible Poisson's ratio. For 
the values of anisotropic material properties chosen, the effect of 
dewetting is to decrease stress levels throughout the propellant grain. 

For the high modulus value of the case used in this calculation, all 
stresses within, the propellant are compressive, whether dewetted or not. 

A second program calculates the stresses developed at various points 
within the propellant as a function of the interface between anisotropic 
and isotropic material per degree Fahrenheit cooling below the stress-free 
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temperature of the motor. For this case all stresses are tensile, and again 
a decrease in both tangential and radial stresses occurs throughout the 
propellant as the anisotropic layer grows. Results are given in Figures 
17 and 18. 

Each of these programs also calculates an equivalent stress for use 
in the reaction rate equation. The particular form chosen to describe 
equivalent stress is discussed in appendix III. 
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APPENDIX III 



CALCULATION OF POSITION OF DEWETTING INTERFACE WITH TIME 



An equivalent uniaxial stress function is chosen in agreement with 
Lindsey's elastic dewetting criterion. From equation (6), 



eq 






° t 3 



3 [- 



Lindsey has reported that for composite propellants, the dewetting stress 
in triaxial tension o has been found to be approximately twice the 
dewetting stress in uniaxial tension o^. With this information. 






eq 6 ' 6 “ 2 

For plane strain conditions: 



o. + o 
6 r 



, 1/2 



°e - °r ' 



For the cases under consideration, > o , thus the absolute value 

0 r 



signs may be removed to yield o 
used during computation. 



eq 



.97cr - .47c . This relation was 

6 r 



Parameters chosen for the reaction rate equation are those obtained 
12 

by Graham and Robinson for rupture of filled PBAN propellant: 

A = -18.1 B - 10104 C « 0.055 



While these values cannot be assumed to be typical of dewetting in PBAN 
propellant or applicable to the entire temperature range, they will be 
satisfactory for sample calculations. 
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For the internal pressurization of a cylinder in steel casing, 
stresses obtained from the elastic solution produce negative values for 
the equivalent stress throughout the cylinder, indicating that dewetting 
is entirely suppressed. For the purposes of calculation, a much more 
flexible case was assumed, with a ratio of moduli 

^propella nt = 0>(J1 
£ 

case 

The time to dewetting at each radial position is calculated with the above 
choice of constants and the geometry described in appendix II. Results 
are given in Figure 19. 

For the thermal cool-down problem, a stress-free temperature of the 
encased propellant was required. This was chosen to be 172°F. A rate of 
cooling of 0.1°F per minute was chosed for calculation purposes. Figures 
20 and 21 show the calculated time to dewetting at each radial position 
as the cylinder and steel case are cooled isothermally . An interesting 
result in this case is that the entire cylinder has dewetted by the time 
the temperature has dropped thirty degrees from the stress-free temperature. 
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